Bogolyubov approximation for diagonal model 
of an interacting Bose gas 



M. Corgini a , D.P. Sankovich b '* 

^Departamento de Matemdticas, Universidad de La Serena, 
Cisternas 1200, La Serena, Chile 

b V.A. Steklov Mathematical Institute, Gubkin str. 8, 119991 Moscow, Russia 



Abstract 

We study, using the Bogolyubov approximation, the thermodynamic behaviour of 
a superstable Bose system whose energy operator in the second-quantized form 
contains a nonlinear expression in the occupation numbers operators. We prove 
that for all values of the chemical potential satisfying fj, > A(0), where A(0) < is 
the lowest energy value, the system undergoes Bose-Einstein condensation. 
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1. Exactly solvable models of strongly interacting bosons could be helpful in 
understanding the nature of Bose-Einstein condensation (BEC) and superflu- 
idity in interacting Bose gases. In this paper we study the thermodynamic be- 
havior of some theoretically relevant system of Bose particles. The considered 
model has rather simplified character. Researches related to similar models 
are connected with the attempt to consider the effect of excitation-excitation 
coupling in Bogolyubov theory. 

The aim of this paper is to consider a model of Bose gas with diagonal interac- 
tion which shows an independent on temperature BEC. This kind of BEC was 
theoretically discovered by Bru and Zagrebnov for some specific Bose system 
with diagonal interactions [Tf2] . Here we study a nonlinear modification of the 
Bru-Zagrebnov model. Furthermore we want to clear up the influence of high- 
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order perturbations in the non-ideal Bose gas interaction on the magnitude of 
the Bose condensate fraction. 

The plan of the paper is as follows. Firstly, we describe the basic mathematical 
notions associated to this kind of systems. Secondly, we discuss the thermo- 
dynamic behavior of a superstable Bose system whose Hamiltonian contains 
a nonlinear term in the number operators. The nonlinearity can be physically 
understood as the simultaneous creation of fc-single one mode particles after 
the disappearance of an equivalent amount of single particles associated to the 
same mode. For k = 2 a variation of this kind of models has been extensively 
studied in refs. [Tf2"f5P] by using different mathematical techniques. (In con- 
trast to the Bru-Zagrebnov model [2J, the Hamiltonian we consider contains a 
specific nonlinear mean-field term (see also [5])). In our case use will be made 
of the so-called Bogolyubov approximation [6]. In our proof a significative 
role plays the fact that the involved Hamiltonian is written in terms of the 
occupation numbers operators. It enables us to give a simplified proof of ther- 
modynamic equivalence of the limit grand-canonical pressures corresponding 
to the energy operator of the system and its respective approximative Hamil- 
tonian. 

We have to point out that although this modified approach leads to a straight- 
forward proof of BEC based on the commutativity of the involved operators, 
it does not apply for proving spontaneous breakdown of symmetry. 

2. Let A be the operator Laplacian. The one-particle free Hamiltonian corre- 
sponds to the operator = — A/2 defined on a dense subset of the Hilbert 
space H l = L 2 (Ai), being Ai = [-1/2, l/2] d C W 1 a cubic box of boundary dAi 
and volume Vi = l d . We assume periodic boundary conditions under which D\ t 
becomes a self-adjoint operator. In this case, as in ref. [3]T7] . we shall consider 
a negative isolated lowest value of energy. This can be done by adding a term 
to the original Hamiltonian which shifts the energy level of the zero mode of 
the kinetic energy operator downward by a negative amount, creating in this 
way a spectral gap. 

We shall study model of interacting Bose particles whose Hamiltonian is given 

as 

^ = ^° + ^t(4« + ^ l t E H) k 4 + v l9 (^-), (i) 

v l v l P eA*\{0} V Vl ) 

where integer k > 2 and g(x) is a suitable real valued function satisfying 
general conditions ensuring the superstability of the model. In particular, it 
suffices to assume that g(x) is a continuously differentiable function on [0, oo) 
satisfying g(0) = and g (x) — > oo as x — > oo [8]. The sum in (JTJ) runs over 
the set Ai = {p e R d : p a = 2Tm a /l,n a = 0, ±1, ±2, . . . , a = 1,2, ...,d}, 
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the Bose operators of creation and annihilation of particles defined 
on the Fock space JF B satisfying the usual Bose commutation rules: [a p , at] = 

5p : q : being S P:q the Kronecker's delta. Denote by 
the occupation-number operator in the mode p G A*. Let Xi(p) = p 2 /2 for 
p ^ and Aj(0) = A(0) < 0. In this case Hf = J2 P eA* Xi(p)n p , 70,71 > 0. 

N — J2 P eA* n p is the total number operator and N' = Y, P eA*\{o} n v * s the 
number operator with exclusion of n§. Note that the boson Fock space JF B is 
isomorphic to the tensor product ® p ^a*^ where is the boson Fock space 
constructed on the one-dimensional Hilbert space H p = {•ye ipx }. 



Let 



be the grand-canonical pressure corresponding to Hi, where p represents the 
chemical potential and (5 = 6~ l is the inverse temperature. 

If Hi(p) = Hi — pN, the equilibrium Gibbs state (— is defined as 



(A) Hl0i) = 1*^,?* exp ( -(3Hi(p) ) 



for any operator A acting on the symmetric bosonic Fock space. The density 
of particles at infinite volume is defined as 

/N\ 

p(p) = Jim pi(p) = Jim ( — ) 

Vj-oo V t ^oo \ Vl / H,(a) 



The density of particles of energy A(0), corresponding to the mode p = 0, is 
given as 

Po{P,fJ>) = Jim Po,i{(3,p) = Jim (^) 

V ; -»oo Vj-oo \Vi / h,(u) 



In this paper we want to prove the occurrence of macroscopic occupation of the 
ground state, i.e. we attempt to find values of p and (3 for which the inequality 
Po(P,p) > holds. Indeed, we prove for the studied model that macroscopic 
occupation independent on temperature takes place, i.e. the density of the 
condensate not depends on j3. 

3. In the framework of the Bogolyubov approximation we must replace the 
operators a\ and a$ in any operator A expressed in the normal form in the 
Hamiltonian Hi(p) with the complex numbers \fVic and \fVic. We get the 
following approximative Hamiltonian 
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P eA*\{0} 



7i s.ukk 



(Xi{p) - n)n p + TFfc=r(oJ) a p 
v i 



+ V W 



+(A(0)- A i)^|c| 2 + 7o ^|c 



2k 



CG C. 



Noting that (a P ) k (a p ) k = n p (n p — 1) • • - (n p — k + 1) for all p G A*, it is not 
hard to see that both Hamiltonians ([1]) and ([2]) are diagonal with respect to 
the occupation-number operators. 

In this case we have the following result [9]. 

Theorem 1 Tr^ , w ~b e -/Jff," ppr (c/0 < Tr , f Be" mw . 

4. Let /v; (x) be the real valued function defined as 



fv,(x) = (ji- A(0))x - 7 x (x - — 



We define the functions 

f L (x) = (fi-X(0))x~ lo x k 



and 



$(x) = <ji-\(0))x-ya[x 



k-l 

V, 



For a; > {k — 1)/Vj the functions f L (x) and fv,(x) are concave with unique 
maxima attained at x L and , respectively. Evidently, 



lim x 



X 



' fc-A(O))- 

7o& 



i 

fc-i 



where (a)+ = max(0, a). 

Lemma 2 For VJ sufficiently large such that x L > (k — 1)/VJ, i/ie following 
inequalities 

f h (x L )< sup /vj(x) < /g(a#,), 

[^p°°) 

. . . A; — 1 (2(Jfc - 1))! 

sup / Vj a; < ^_A(0) +— — + 7o- 
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hold. 



PROOF. For V\ sufficiently large and x G [(k — 1)/VJ, oo) we have the obvi- 
ous inequalities f L (x) < fv t (x) < fv,( x )' which for x L > {k — 1)/VJ imply the 
first statement of the proposition. On the other hand, 



f Vl {x) < (// - A(0))+x + 7 x (x + 



X + 



k- r 



v. 



for x E [0, (k - l)/Vi). Then, 

sup / Vi (x)<(//-A(0)) + +7o 



□ 



Lemma 3 T/ie following inequality takes place: 



lim sup f Vl (x) < f (x 



V,^oo 



[0,oo) 



fc-A(O))- 
^7o 



(A*-A(0)). 



fe-1 



PROOF. 



l/?(*g) "/V)l = /vi K + 



V5 



/ L (x L ) = (/x - A(0))+ 



fe-1 

~1T : 



i.e. limy^oo /^(xK) = / L (x L ). Then, using this result and the first inequality 
in lemma [2j we get 



L/„,L\ 



lim sup f Vl (x) = f (x 

Vi—>oo ,k-l \ 



fc-A(O))- 
&7o 



Gu-A(0)). 



fe-1 



The second inequality of lemma [2] yields 



lim sup fv t (x) < 0. 



[o, V) 



From these results, noting that 



sup f Vl (x) < max{ sup fv t (x), sup f Vl (x)}, 

[0,oo) [„,*=!) [*=!,«,) 
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and passing to the limit Vi — > oo, the proof of the lemma follows. 
□ 

Let p — n /Vi. p is a densely defined on JF^ positive self-adjoint operator with 
spectrum cr(p ) C [0, oo) and spectral family {Ep (x) : x G o"(po)} [10J . Let 
: [0, oo) — ► R be a continuous function such that sup^o^ h(x) exists. 
Assume also that h(x) exp(pVifv l ) is a bounded and continuous function in 
[0,oo). Let A, B be self-adjoint bounded operators in the Hilbert space U 
endowed with inner product denoted as (■,■)• Then we say that A is smaller 
than B and write A < B if (0, A<f>) < ((f), B<f) for every G W. 

Lemma 4 (HPo)) HiQj.) < sup /i(x). 

a;G[0,oo) 



PROOF. The spectral theorem for self-adjoint operators [TU] enables us to 
write 

h{po) exp(pVif Vl (po)) = J h(x)exp(pVif Vl (x))dE po (x). 

o-(po) 

Although po is an unbounded operator, the above representation is valid for 
any G J-~q since h(x) exp^Vifv^x)) is a bounded continuous function. This 
implies that, 




in the sense of operators for all G Jq . Noting that 
exp(-f3Hi(p)) = exp(^/v,(p ))exp(-/3^'(/i)) 

and since exp(/?Vz/vi(po)), exp(— f3H{(p)) are positively definite operators com- 
muting with each other, it is not hard to see that 

h(po) exp(f3Vif Vl (p )) exp(-/3i^'(//)) 

< sup h(x)exp((3V l f Vl (p ))exp(-(3H' l (p)). 

x<={0,oo) 

Then, taking into account that Hi(p) defines a superstable system and since 
JF B is isomorphic £g> pe A ; * ^ we get 
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Tr^B h(p ) exp(/3VJ/y ; (p )) exp(-/3fl,'(/z)) 

< sup h(x)Tr^ B exp((3V l f Vl (p ))exp(-(3H' l (fi)). 

x<={0,oo) 

Hence, we conclude that 

(Hpo))hu) < sup h(x). 

xe[o,oc) 

□ 

Theorem 5 The systems with Hamiltonians Hi(p) and iJf ppr (c, fi) are ther- 
modynamically equivalent in the Wentzel sense /IT) /, i.e. 

lim sup ^ ppr (/3,c,/i) =p(fi,fi), 

v i^°° \ c \:ceC 



where 



P f ppr (Ac,/i) = ^lnTr^ An{o} ^exp[-/5(^ ppr (c )/U ))] 



andp(f3,fi) = lim Vl -+o C Pi(/3, p). 



PROOF. Note that this is an equivalence of the Wentzel type but considering 
infimum for the real variable |c| instead for c G C. Applying the Bogolyubov's 
convexity inequalities for pressures and theorem [1] we obtain that 

< Pl ((3,p:) -pTO?.^) < Unr\c^) -Hfo))^), 

v\ 



where 



hHr T {c^)-H^)) HlM = (A(0)-/i)|c| 2 + 7o | 

Vi 



C 



+ / - A(0))p - loPo ■■■(po- 



Blip) 



Hence, we obtain 

< Pl ({3,ti)-p? ppi (P,c,ti) < (A(0) - ^\c\ 2 + l0 \c\ 2k + (f Vl (p )} HlM . (3) 
We have to distinguish two cases. 
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(i) Case ji > A(0). In this situation, inf| c |. ce c (A(0) — fi) \c\ 2 + 7 |c 
tained for \cq\ given by 



2k 



is ob- 



l c o| 



>-A(0) \ 



This leads to 



inf 

|c|:ceC 



(A(0) - //) \c\ 2 + 70 \c\ 2k 



&7o 



k 

f h {x h ). 



From the above result, using lemma H] with = fvi( x ) and considering the 
infimum for c G C in the right-hand side of inequality ([3]), it follows that 



o< Pl ((3,ri- sup pr\^c^)<-(^ 1 ^Y 1 ^-m) k -r^ 

\c\-.cec \ k-f J k 

+ sup f Vl (x). 
[0,00) 

Therefore, passing to the limit V\ — > 00 and applying lemma [3j we get 
P(A p) - ,) im SU P ft Ppr (/^, c, /i) = 0. 

^i^ 00 |cj:cgC 



1 2A- 



is attained 



(ii)Case < A(0). In this case inf| c | :ceC (A(0) — ji) \c\ +7o|c| 

at I Co I = 0. Moreover, since (a\) k aQ is a positive operator, for fi < A(0) we 
obtain 



(/i - A(0))p - 7oPo • ■ • Po 



k - 1 



(fj, - A(0))p - 7o 



(4) fc Qo \ 



< 0. 



Then, from inequality ([3]) we get, 
Pi(P,H) = sup pf ppr (/3,c,/i) 

|c|:ceC 



for any finite VJ. Therefore, in both cases, passing to the limit Vi — > 00, we get 
p(ft P) = J im SU P pf Ppr (/5, c, p). 

|c|:ceC 
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The proof is complete. 

□ 



5. The following theorem holds. 

Theorem 6 For all \x > A(0) the system (QP displays independent on temper- 
ature BEC and the amount of condensate is given by 



PROOF. Using the convexity of p/(/3, /i) and p^ ppr (P, c , //) with respect to 
A(0), we get from a Griffiths theorem [12] 

PoM = -d\(o)Pi(P,fJ>) = -^A(o)Pz appr (/3,c ,/i) = |c | 2 . 
Finally, passing to the thermodynamic limit, we obtain for all /i > A(0), 



6. For the class of Hamiltonians given in equation (PQ), in the framework of 
the so-called Bogolyubov approximation [6|9] . it has been given a simple 
proof of thermodynamic equivalence of the limit grand-canonical pressures 
corresponding to those systems and their respective approximating ones for 
any integer k > 2. Moreover, in contrast to the Bru-Zagrebnov models [TP] , 
we prove that independent on temperature BEC in the sense of macroscopic 
occupation of the ground state holds for any integer k > 2 and any \i > A(0). 
A similar type of BEC is explained entirely by superstability of the model and 
by absence of an interaction between the ground state occupation number 
operators and the nonzero modes ones. 
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